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ABSTRACT 

We investigate expected sub-mm emission and dust properti es of the Lyman Brea k Galaxies (LBGs) 
in the Hubble Deep Field (HDF). The SCUBA Deep Survey ( [Hughes et al., 1998Q provides an upper 
limit of the 850/im flux densities of the LBGs. With this constraint, we argue that a typical ultra-violet 
to far-infrared spectral shape of the high-redshift LBGs is rather close to a template spectrum of low- 
reddening starburst galaxies in the local universe but different from that of heavily dust-enshrouded 
ultra-luminous far- infrared galaxies like Arp220. We also evaluate the lower-limit temperature of dust 
in LBGs assuming single- and two-component modified black-body spectrum. To estimate the total 
amount of energy re-emitted in FIR wavelength, we take two different approaches; model fitting of the 
UV spectra of LBGs and an empirical method using the relationship between UV spectral shape and 
UV/FIR flux ratio observed for local starburst galaxies. Both methods give lower-limit temperature of 
the LBGs as ~ 40K, which is higher than the typical dust temperature of local optical- and infrared- 
selected galaxies. This result is also supported by the comparison of the expected sub-mm flux of the 
LBGs with the cosmic FIR background radiation. The high dust temperature may indicate the effective 
massive-star formation or different dust properties in the high-redshift LBGs. 

Subject headings: galaxies: starburst — galaxies: ISM — galaxies: formation 



1. INTRODUCTION 

A number of high-redshift star-forming galaxies have 
been identified through the Lyman break surveys (Stei- 
del et al. 1996a, b; Lowenthal et al. 1997). Since the 
number density of these Lyman Break Galaxies (LBG) is 
comparable to or somewhat higher than that of the local 
~ L* galaxies, they are thought to be natural candidates 
of the progenitors of nearby normal galaxies (Lowenthal 
et al. 1997). From the observed UV continuum emis- 
sion, a typical star-formation rate (SFR) of LBGs is es- 
timated to be moderate, ~ 10 M Q yr _1 (Ho=50 km s _1 
Mpc -1 and go=0-5). However, there are increasing evi- 
dences of the presence of significant amount of extin ction 



in the UV continuum by interstellar du st in LBGs (|Pct 
tini ct al 1998fc |5awicki fc Yee 1998[ ). If the inferred 



UV extinction is corrected, their SFR would be more than 
£ 100M© yr" 1 (Dickinson 1998), which may imply that 
high-redshift LBGs are in their major star-formation pe- 
riod. 

On the other hand, deep surveys in sub-mm wave- 
length with Sub-millimeter Common-User Bolometer Ar- 
ray (SCUBA; Holland et al. 1998) revealed a population of 
FIR-luminous galaxies possibly at high redshift and have a 
large SFR, £ 100 M yr" 1 (Smail et al 1997; Hughes et al. 
1998; Barger et al. 1998; Lilly et al. 1999). Although the 
optical identification of the SCUBA sources has not been 
well established in some cases, they are also candidates of 



high-redshift forming galaxies. 

It is then natural to ask what the relationship between 
these two populations is, or, how the high-redshift LBGs 
appear in sub-mm wavelength. In the deepest sub-mm 
imaging so far at the Hubble Deep Field (HDF; Williams 
et al. 1996) by Hughes et al. (1998), none of the optically- 
selected LBGs (Steidel et al. 1996b; Lowenthal et al. 1997) 
was detected. The upper-limit flux value at 850 //m, how- 
ever, can be used to put some important constraints for the 
spectral energy distribution (SED) of the LBGs. In this 
Letter, we investigate sub-mm emission properties of the 
LBGs in HDF to constrain their physical properties such 
as dust temperature by using the results of the SCUBA 
observations by Hughes et al. (1998). 

2. SPECTRAL SHAPE OF THE LYMAN BREAK GALAXIES 

Our sample of LBGs is the same as that used in Sawicki 
and Yee(1998). They are the 17 galaxies at 3 with firm 
spectroscopic redshifts and with photometric data from 
optical to near-infrared (NIR) wavelength, taken from the 
41 high-redshift LBG candidates originally provided by 
Steidel et al. (1996b) and Lowenthal et al. (1997). Red- 
shifts and magnitudes are listed in Table 1 of Sawicki and 
Yee (1998). 

Hughes et al. (1998) detected five sources at 850 /^m 
with SCUBA in HDF. We found none of the optically- 
selected LBGs is detected within the error circles of the 
SCUBA sources 3 . Hughes et al. evaluated their detec- 
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tion limit ~ 2 mJy (5a). Hereafter we quote this value as 
the upper-limit flux density at 850 fim for all the LBGs al- 
though the true detection limit may change across the field 
slightly due to the position-dependent confusion noises and 
detector sensitivities. 

The upper limit of sub-mm flux density gives a con- 
straint on the spectral shape from the rest-frame UV to 
FIR wavelength of the LBGs. We compare the observed 
optical-NIR flux densities and the upper limit of sub-mm 
flux density of LBGs with the three template spectra of lo- 
cal starburst galaxies. For the template spectra, we adopt 
the SED of Arp220 (Klaas et al. 1997), as an example 
of heavily dust enshrouded starburst galaxies, and the av- 
eraged SED of high-reddening starburst galaxies (SBH) 
and low-reddening starburst galaxies (SBL) compiled by 
Schmitt et al. (1997). These templates are redshifted and 
then fitted to the optical-NIR photometric data. Figure [l]a 
and |l|b show the best and the worst case in the results of 
our least-square fitting procedure. In many cases, the tem- 
plate of SBL gives a better fit but those of Arp220 and SBH 
are always steeper than the spectra of LBGs (Figure [TJa). 
For several objects, spectra of LBGs are even flatter than 
the SBL template (Figure |l|b). Reduced \ 2 values (de- 
gree of freedom is four) ranges 300-550 for Arp220, 50- 
250 for SBH, and 1.5-30 for SBL. The absorption at the 
rest-frame UV wavelength in the LBGs thus seems much 
smaller than in Arp220- or SBH- like galaxies. If we assume 
the templates of Arp220 and SBH (fitted at optical-NIR 
wavelength), the expected flux density at 850/im exceeds 
2 mJy for 11 and 9 of the 17 objects, respectively. If the 
SBL template is applied, none of the objects has a 850/im 
flux density larger than 2 mJy, which is consistent with 
the results of the SCUBA observation. Thus, not only the 
rest-frame UV-optical SEDs but also the FIR-UV SEDs of 
LBGs are better represented by the SBL template. 

On the other hand, SCUBA-selected bright sub-mm 
sources show more FIR excess. For example, Barger et 
al. (1998) detected a sub-mm source with S^sso/im) = 4.6 
mJy by SCUBA in the Lockman-Hole field and found a 
plausible optical counterpart with Kab = 21.8. They ar- 
gued that only the SED of ultra- luminous infrared galaxies 
like Arp220 can match this flux ratio. The SEDs of the 
LBGs are thus much shallower than those of the SCUBA 
sources detected above a few mJy level at 850 /im. The 
observed star-forming regions of LBGs are not likely to be 
heavily enshrouded by dust compared with the SCUBA- 
selected bright sub-mm sources or the local ultra-luminous 
infrared galaxies. 

3. DUST PROPERTIES OF LYMAN BREAK GALAXIES 

3.1. Model Fitting of the UV Spectra 

We first estimate the amount of dust extinction of LBGs 
by fitting the model UV/optical spectra (at 900-10000 
A; Salpeter IMF is assumed ) with an extinction to the 
observed SEDs following the similar manner as in Saw- 
icki and Yee (1998) but using the different evolutionary- 
synthesis model developed by Kodama & Arimoto (1997). 
We assumed the reddening law derived by Calzetti for 
the local starburst galaxies 4 . The estimated extinction, 
E(B — V), ranges 0.3-0.5 mag with ages of ~ 0.03 Gyr, 

LBGs discussed here. 

4 The formula is presented in Sawicki and Yee (1998) 



which agrees with the results in Sawicki & Yee (1998). For 
13 of the 17 objects, acceptable fittings are obtained with 
finite amounts of extinction and we will discuss the 13 ob- 
jects in this section. The 850/im flux densities of LBGs 
are then calculated, with the assumption that all the en- 
ergy absorbed in UV/optical wavelength is re-emitted by 
the modified black-body radiation with single tempera- 
ture. The observed 850 /im flux density, S^ssrj/jm) ; is given 

by 

(1 + z)F a i, s g 
Jo vPB v (T)dv 

where F a ha is absorbed flux, T is a dust temperature, 
(3 is a power of the emissivity law, B V (T) is the Planck 
function, and v x = (1 + z)i/g 5QlJ , m . 

Figure |^ shows the dependence of the estimated flux 
density on the dust temperature and emissivity for a typ- 
ical case. We assume (3 = 1 in the following discussion, in 
order to compare with the properties of the local galaxies 
which are obtained with (3=1. For the given (3 value, 
the upper limit of the flux density at 850 /im (2 mJy) im- 
poses the lower limit of the dust temperature, Tu m . In 
Figure |^ (the fourth panel from the top), the distribu- 
tion of the obtained Tu m is shown along with those of 
the dust temperature of local galaxies from the samples of 
the optically-selected CfA galaxies classified as E-SO/a and 
Sb-Sbc (Sauvage and Thuan 1994) and the I R AS -selected 
infrared-bright galaxies (Young et al. 1989). These dust 
temperature of the local galaxies are measured by fit- 
ting the IRAS 60 and 100 /im flux densities with single- 
temperature modified-Planck function with f3 = 1. The 
distribution of Tu m of the LBGs has a peak at around 45 
K and significantly higher than those of the local galax- 
ies. The median value of the temperature lower limit 
is 43 K while those of CfA and infrared-bright samples 
are T med = 33 K (E/S0), 30 K (Sb-Sbc), and 36 K (in- 
frared bright). According to the Peto-Prentice General- 
ized Wilcoxon Test, the temperature distribution of LBGs 
is different from those of the local galaxies with a confi- 
dence level larger than 99.9999 % . 

It is possible that the high dust temperature of the LBGs 
is introduced by our assumption of single dust temper- 
ature. Most of the nearby star-forming galaxies have a 
hot dust component (T ~ 100 — 150 K) in addition to a 
cold one (T ~ 30 - 50 K) (Eales, Wynn- Williams, and 
Duncan 1989; Klaas et al. 1997). IRAS 60 and 100 /im 
flux densities of the local galaxies may be dominated by 
only the cold component. So as to estimate the energy 
drained to the hot dust, we fitted the SEDs of the three 
template spectra used in the previous section by the two 
component models. The results are L(cold) / L(hot) = 1.9 
for the SBL, L(cold) / L(hot) = 3.2 for the SBH, and 
L(cold) I L(hot) = 9.8 for the Arp220 template (The result 
for Arp220 is consistent with that of Klaas et al.). If we 
consider the contribution by the possible hot component 
in calculating the FIR luminosity of the LBGs, adopting 
the ratio L(cold) / L(hot) = 1.9 which gives the lowest lim- 
iting temperature, the peak of Tu m of the LBGs is shifted 
to ~ 37 K. According to the Peto-Prentice Generalized 
Wilcoxon Test, they are still different higher than 99.9999 
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% level. True Tu m must be between these two extreme 
cases, and it may be Ti im ~ 40 K. In any case, the dust 
temperature of the LBGs are likely to be still higher than 
those of the nearby normal galaxies. 

The difference of local galaxies and high-redshift LBGs 
can be in dust extinction properties instead of dust tem- 
perature. The reddening law in the LBGs is not under- 
stood very well (e.g., Pettini et al. 1998). If we use the 
SMC-like reddening law instead of the Calzetti law, the 
total amount of absorbed UV luminosity becomes smaller 
and the expected sub-mm flux density becomes several 
times smaller. In such a case, only a weak constraint for 
dust temperature is obtained. 

3.2. FIR emission Estimated from UV spectral Slope 

UV spectra of the local starburst galaxies are well ap- 
proximated by a power-law of index p (1250 ~ 2200 A), 
j a = A\ p (ICalzctti Kinney fc Storchi-Bcrgmann 1994, 



Meurer et al., 1995| ), and there is an empirical rela tion 
between p and FpinJ 'Fjjv values ( Meurer et al., 1997| ). In 
the range of p ^> — 1.8, approximately, 



i? ee log 



Ffir 
Fuv 



0.65p+ 1.8 



(2) 



where F FIR = 1.26 xlO -11 



2.58 



(^H + (- 



Jy 



and F uv = A c / Ac (A c = 2320 A). Hence, 
Ffir « AX/+ 1 x 10*. (3) 



If we assume the modified black-body radiation with the 
power of emissivity law (3 = 1.0, the dust temperature T 



implies the flux ratio k(T) = S^ioo^m) 
we obtain, 



IS v (60 fi 



1.26 x 10- n (2.58 + k{T)) Vergs" 1 cm" 2 



m) • 



FIR 
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The sub-mm flux is predicted by the spectrum from eq. 
Again, the upper limit of S v Km(850/im) = 2 mJ y S ives T Um 
values. The bottom panel of Figure [3] shows a distribution 
of Tu m obtained by this empirical method. It peaks at 
~ 38K, which is similar to or somewhat lower than the 
value calculated from the model UV extinction. It is not 
surprising that the empirical method gives wider Tu m dis- 
tribution since there is fairly large scatter in the empirical 
relation (eqj|) itself (Meurer et al 1997). 

4. CONTRIBUTION TO THE CFIRB 

The average dust temperature of LBGs can also be ex- 
amined with the observed intensity of cosmic far-infrared 
background (CFIRB) at 850/im, although the field of view 
of HDF (4.7 arcmin 2 ) is very small and statistical un- 
certainty may not be negligible; a typical fluctuation of 
CFIRB at this scale is still unknown. The total 850/im 
flux in HDF evaluated from the average intensity of the 
CFIRB radiation, vl v = 3.5 x lO" 10 W m~ 2 sr" 1 (Fixsen 
et al. 1997). The five SCUBA sources contribute about 



45 % of this flux. The contribution of the LBGs in HDF 
should not significantly exceed the remaining 55% of the 
CFIRB flux. 

In Figure [|, we plot the contribution of the LBGs stud- 
ied in the previous section to the CFIRB with respect to 
the assumed dust temperature. Three curves are those 
given by the model calculations with single and double 
temperature components and the prediction using the em- 
pirical relation. Depending on temperature, LBGs con- 
tribute a significant fraction (^> 30%). With the constraint 
not to exceed the CFIRB flux, we obtain fairly high lim- 
iting temperature, Tu m ~ 40 K, which is consistent with 
the limit obtained by using the SCUBA upper-limit flux 
density. 

5. DISCUSSIONS 

We showed that the LBGs in HDF have UV-to-FIR 
SEDs different from those of the ultra-luminous FIR galax- 
ies in the local universe and of the SCUBA-selected bright 
sub-mm sources. It may imply that the birght sub-mm 
sources detected above a few mJy at 850 /im are rather 
extremely dust-rich objectes among the high-redshift star- 
forming galaxies selected by the Lyman break technique 
although it does not have to mean that the entire parent 
populations of sub-mm sources and LBGs are fundamen- 
tally different. 

We also claimed that dust temperature in the LBGs 
should be higher than ~ 40K, if we assume the reddening 
law appropriate for local starburst galaxies. What causes 
the relatively high dust temperature in the LBGs? High 
dust temperature may be results of some intrinsic proper- 
ties of star formation in high-redshift galaxies. There are 
some possible explanations for effective dust heating. For 
example, a flatter initial mass function (IMF) which pro- 
duces a larger number of massive stars for a given mass 
of dust may provide the effective dust heating. High star- 
formation efficiency (a higher SFR for a given mass of dust 
regardless of IMF) brings the similar situation. The dif- 
ference may also be in properties of dust grains. For the 
case of low metal abundance (possibly true in young galax- 
ies), dust may be dominated by small-size grains which are 
easily heated. 

Finally, we discuss the detectability of LBGs in sub-mm 
observations. Assuming an appropriate spectral shape like 
SBL or dust temperature of 40 — 50X, we predict that 
sub-mm flux at 850/im of the LBGs in HDF is typically 
0.1-1 mJy. So within the reasonable amount of observa- 
tional time, only the brightest LBGs could be detected 
with SCUBA. In the coming decade, however, new sub-mm 
observational facilities such as Large Millimeter and Sub- 
millimeter Array (NAOJ) or MilliMeter Array (NRAO), 
etc. are planned. Our results presented in this letter pro- 
vides rather optimistic future prospects; only a few hour 
exposure time will be needed to detect typical LBGs with 
such new facilities and then we will learn about true star- 
formation properties of typical high-redshift galaxies. 

We thank N.Arimoto and T. Kodama for kindly provid- 
ing their evolutionary synthesis code. 
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Fig. 1. — SED of the LBGs is fitted by the templates of local starburst galaxies; high-reddening starburst (dashed line) and low-reddening 
starburst (solid line) compiled by Schmitt et al. (1997) and Arp220 (Klaas et al. 1997; dash-dotted line). Examples of the best and the worst 
fittings arc shown in panel (a) and (b), respectively. The [/-and B-band data are not used in the fitting. 
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Dust Temperature [K] 



Fig. 2. — An example of the dependence of the estimated flux density at 850 (im on the dust temperature and the power of dust emissivity 
law p. 
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Fig. 3. — The bottom two panels show the distributions of the lower-limit dust temperature of LBGs evaluated by modeling the UV 
extinction (single dust temperature assumed) and by using the empirical method (see text). The top three panels show the dust temperature 
of the local optically- and infrared-selected galaxies in literatures. 
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Fig. 4. — The expected contribution of LBGs to the CFIRB at 850^tm is plotted in function of the assumed dust temperature. The cases of 
the model extinction with single- and two-component dust emission (solid and dashed line, respectively) and the case of the empirical method 
(dash-dotted line) are shown. 



